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Abstract-Equations  by  Soulsby  (Dynamics  of  Marine  Sands, 
Thomas  Teleford:  London,  1997)  and  Whitehouse  (Scour  at 
Marine  Structures,  Thomas  Teleford:  London,  1998)  have  been 
applied  to  predicting  mine  burial  by  scour  in  tidal  estuaries.  The 
main  equation,  which  has  the  functional  form  of  [1-  exp(-t/T)], 
computes  the  depth  of  the  scour  pit  with  time  under  steady  state 
conditions.  This  theory  may  be  applied  to  changing  conditions 
by  using  the  RMS  values  of  the  frictional  stress  at  the  bed  and 
assuming  a  quasi-steady  state  of  the  RMS  values  over  a  small 
time  period.  Lrom  1999  -  2002,  the  Naval  Research  Laboratory 
conducted  scour  burial  experiments  using  instrumented  mines 
that  measure  mine  motion  (heading,  roll,  and  pitch)  and  percent 
burial  (surface  area  covered  with  sediment).  Using 
oceanographic  and  sediment  data  obtained  during  these 
experiments,  this  study  examines  how  well  the  predictions  match 
the  mine  burial  measurements. 

I.  INTRODUCTION 

The  Navy  has  recently  shown  increased  interest  in 
predicting  mine  burial,  both  because  of  mining  activities  in 
the  Persian  Gulf  and  because  the  state  of  prediction  has 
arguably  been  poor  [1],  For  sandy  environments,  mine  burial 
often  occurs  from  scour  -  the  erosion  of  sediment  around  an 
obstacle  due  to  intensified  shear  stresses  and  vortices  caused 
by  the  object’s  presence.  (Other  burial  processes  can  include 
sand  ridge  migration,  liquefaction,  and  global  sediment 
transport;  impact  burial  usually  is  insignificant  for  sandy 
bottoms.)  The  most  commonly  used  scour  mine  burial 
models,  until  recently,  were  based  on  equations  presented  in  a 
1963  report  by  Carsten  and  Martens  [2].  Their  measured  rates 
of  burial  are  based  on  scaled  laboratory  experiments,  which 
may  not  represent  in  situ  conditions  where  the  mines  are 
much  larger  than  mean  grain  size  or  rippled  bed  forms. 

A  newer  and  more  promising  scour  burial  model  has 
recently  been  identified  by  Friedrichs  [3].  The  model  comes 
from  a  set  of  equations  developed  by  HR  Wallingford  Ltd.  (a 
United  Kingdom  civil  engineering  firm)  and  presented  in 
books  by  Soulsby  [4]  and  Whitehouse  [5].  It  was  developed 
to  predict  scour  in  steady  state  conditions,  but  has  been 
adapted  to  changing  environments.  This  adaptation  is  made 
by  assuming  a  quasi-steady  state  of  root-mean-squared 
bottom  orbital  velocities  and/or  bottom  currents  over  a  small 
time  period  and  evaluating  model  output  in  a  time-stepped 
fashion. 

Tests  of  this  model  have  been  performed  using  data 
collected  from  recent  in  situ  experiments  with  an 
instrumented  mine  capable  of  recording  burial  over  a  period 
of  a  few  months.  Data  from  one  of  these  deployments  was 
used  by  Friedrichs  [3]  to  call  attention  to  the  potential 
usefulness  of  the  model.  Comparison  with  data  from  two 
other  deployments  has  shown  similarly  positive  results.  This 


paper  will  present  the  theory  for  this  model,  followed  by 
comparison  of  the  model  output  with  experimental  data  as 
applied  to  burial  of  a  cylindrical  free -body  in  shallow  sandy 
oceanic  environments. 


II.  THEORY 

The  depth  of  the  scour  pit,  S(t),  for  steady  state  flow  after 
an  amount  of  time,  t,  has  been  determined  [6]  empirically  to 
be 


S(t)=Saa{l-exp{t/T)p),  (1) 

where  .S'x,  is  the  scour  pit  depth  after  an  infinite  amount  of 
time,  T  governs  the  rate  of  scour  pit  growth,  and  p  is 
determined  by  the  geometry  of  the  scouring  body.  For  a  5:1 
cylinder,  p  =  0.6  [7].  The  other  two  variables  ( S.^ ,  T)  are 
calculated  from  oceanographic  conditions  and  physical 
properties  of  the  sediment.  The  next  two  subsections  describe 
the  sequence  of  calculations  required  to  obtain  SrJ  and  T  at 
each  time  step,  using  the  statistics  of  the  oceanographic 
conditions  for  that  time  period. 


A.  Calculating 

1 )  Basic  Relation 

Whitehouse  [7]  gives  the  value  of  Srf.  as  a  relationship 
between  the  velocity  at  the  sediment  bed,  U,  relative  to  the 
velocity  required  for  the  sediment  to  become  mobile,  Ucn  (the 
critical  velocity)  in  the  following  manner 


S  oo  —  1 


0, 

l.l5D0(2U-l.5Ucr)/Uc 

IA5D0, 


U  <0.75f/cr 
0.75 Ucr  <U  <  l.25Ucr 
U  >l.25Ucr 

(2) 


where  D0  is  the  diameter  of  the  mine. 

It  is  more  convenient,  however,  to  rewrite  (2)  so  that  the 

9 

magnitude  of  the  stress  at  the  sediment  bed,  r  =  pwCrU  ,  is 

compared  to  the  critical  stress,  zcr  =  pwCrU cr ,  required  for 
the  sediment  to  become  mobile.  (In  these  relations,  pw  is  the 
density  of  water  at  the  sediment  interface  and  C,  is  a 
dimensionless  mechanical  resistance  constant  related  to  the 
roughness  of  the  seafloor  and  strength  of  the  flow.)  The 
reasons  for  this  change  are  that  the  critical  stress  has  been 
deduced  from  laboratory  experiments,  the  stress  at  the  seabed 
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can  be  calculated  from  currents  and  oceanographic 
conditions,  and  the  critical  velocity  has  to  be  obtained  from 
knowledge  of  Cr,  which  is  difficult  to  assess  for  U  =  Ucr. 

The  stresses  are  evaluated  in  terms  of  the  non- 
dimensional  Shields  parameter,  9 , 


The  quantity  zm  is  the  mean  shear  stress  induced  by  waves 
and  currents.  (Equation  (5)  accounts  for  nonlinear 
interactions  between  z  w  and  rc ;  zm  is  in  the  same  vector 

direction  as  T c .)  In  (5),  zc  is  calculated  from  the  assumption 


d  =  Tl((Ps-  Pw)gd5Q>)i  (3) 


zc=pwCDU2,  (9) 


where  ps  is  the  average  bulk  density  of  the  sediment,  g  is 
gravitational  acceleration,  and  d50  is  the  median  diameter  of 
the  sediment  particles.  For  r  =  rcr ,  the  critical  Shield's 
parameter  is  6cr .  It  can  be  seen  from  the  above  relationships 
between  stress  and  velocity  and  stress  and  the  Shield's 

parameter  that  9^ ~ / ' o\!^  =U /U cr  . 

Rewriting  (2)  in  terms  of  the  Shield's  parameters,  one 
obtains 


5 


00 


0,  ^  ^  49ldcr  <a75 

-1.1500(27^7-1.5)  0.75  <  7 0/0cr  <  1-25  , 
1.15D0,  7 6lecr  >L25 


(4) 


The  critical  Shield’s  parameter  has  been  evaluated 
empirically  to  be  [8] 


9cr  — 


0.3 


1  +  1.2D* 


-  0.055(1  -exp(-0.02D*))  (5) 


where 


D*  =  ^50  — 


(6) 


s  =  ps  / pw  ,  and  v  is  the  kinematic  viscosity  of  seawater. 


2)  Obtaining  9 from  Oceanographic  Conditions 
The  last  remaining  unknown  in  (4)  is  9 .  Soulsby  [9] 
provides  the  following  methodology  for  obtaining  total  stress, 
r ,  which  gives  9  by  (3).  The  total  stress  is  obtained  from 
the  stresses  induced  by  wave  motion,  z  w ,  and  by  current, 
rc .  The  following  equation  relates  these  stresses  to  the  total 
stress  at  the  seabed. 


r  = 


+  TX 


,  cost/))  +  ( rw  sin  <fY 


10.5 


(7) 


where  U  is  the  depth  averaged  current  velocity  evaluated 
from 


U(-)=  \{z/032h)V1  U ,  0<z<0.5h 
[l.07(7,  0.5/r  <  z<  h 


(10) 


and  C D  is  the  drag  coefficient  given  by 


CD  = 


0.4 

1  +  ln{\  2h/d5() ) 


2 


(ID 


where  h  is  water  depth.  The  other  stress,  rw ,  is  assumed  to 
have  a  form  identical  to  (9).  Depth  averaged  current  velocity 
is  replaced  with  the  orbital  velocity,  U /, ,  and  the  drag 
coefficient  is  replaced  by  a  wave  friction  factor,  0.5 fw,  so 
that 

zw=0.5pwfwU2,  (12) 

where 

fw=L39-{6UbTw/mi5o)~0-52,  (13) 

and  Tw  is  the  wave  period. 

B.  Obtaining  T 

Once  Soois  evaluated  for  the  current  time  step  from  (2), 
T  is  the  only  remaining  variable  in  (1)  needing  to  be 
evaluated.  An  empirical  relation  for  T  has  been  found  to  be 
[6,  10] 


where 


T  =  T*((s-l)gd350)~U2D2 


T*=A9%. 


(14) 


05) 


where  tf>  is  the  angle  between  the  current  stress  and  wave 
stress  vectors  and 


f  111  —  z  ( 


1  +  1.2 


\3.2 


yTc  +  Tw  j 


(8) 


A  =  0.095  and  B  =  -2.02  for  a  free  cylinder  [7],  The 
co  subscript  on  the  Shield's  parameter  means  that  this 
quantity  is  evaluated  for  when  the  mine  is  absent.  Placing  9 
in  for  9,/;  in  (15)  and  putting  the  result  into  (14)  produces 

T=A9B^s-l)gd350]~l/2D2.  (16) 
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C.  Time  stepping  S(t)  to  calculate  mine  burial  and  re¬ 
exposure 

Now  that  (1)  can  be  evaluated,  it  needs  to  be  re¬ 
expressed  in  a  time-stepped  fashion  to  make  it  applicable  to 
quasi  steady-state  conditions.  The  procedure  developed  by 
Friedrichs  [3,  11]  will  be  used  for  this  purpose.  At  the /h  time 
step,  the  quasi  steady-state  starts  at  t=1!  and  ends  at  t=tj  +  At 
=  tj+i,  where  At  is  the  size  of  the  time  step.  At  t=tj  ,  the 
values  for  S(t),  T  and  are  S(tj),  T(t  j)  and  S^itj). 

Thus,  (1)  becomes 

S(tj+ 1)  =  ^  (tj  )(]-exp(-  ( t  +  At )/T (t  j  ))P ).  (17) 

The  appropriate  value  to  use  for  1  is  the  time  it  would  take  for 
S(t)  =  S(tj)  under  a  steady  state  condition  with  zero  initial 
scour,  T  =  T(tj)  and  .S'x  =  Sm(t  j) . 


where  a  is  the  efficiency  of  the  scour  process  to  expose  the 
mine,  assumed  to  be  0.6  [3]. 

in.  COMPARISON  WITH  EXPERIMENTS 

Three  experiments,  one  at  Scripps  Institute  of 
Oceanography  in  1999  and  two  at  Martha’s  Vineyard  Coastal 
Observatory  (MVCO)  in  2001  and  2002,  have  used  an 
instrumented  mine  (Fig.  1)  to  measure  burial  by  scour  and 
other  post-impact  processes.  Environmental  and  geotechnical 
data  required  to  run  the  scour  model  were  obtained  during 
these  experiments.  Both  sites  had  fine  sand  bottoms  (median 
grain  sizes:  Scripps,  190  microns;  MVCO,  180  microns). 
This  section  presents  results  from  these  measurements  and 
model  predictions.  The  specifics  of  the  equipment  used  and 
the  numerical  methods  employed  will  be  discussed  first, 
followed  by  measurements  and  modeling  results. 


t  =  ~T(tj  )(ln(l  -S(tj )  /  S*  (tj  )))17  P  (18) 

This  condition  calculates  the  process  along  the  appropriate 
portion  of  the  S(t)  curve.  Equation  (18),  however,  assumes 
that  S(tj)  <  S 30 (t  j ) .  (S(t  j)  >  Sro(t  j)  gives  a  complex  time, 

which  is  physically  unmeaningful.)  This  condition  will  be 
false  when  the  scouring  current  is  weakening  because  S(tj) 
was  calculated  at  the  end  of  the  previous  time  step 
while S^itj)  is  applicable  to  the  present  time  step.  In  this 

instance,  S(t j+\)  =  S^it  j)  is  set.  Hence,  (18)  is  rewritten  as 
follows 


S(t j+l)  - 


Sao(tj) 


'  7  (tJ+At)^ 


1-exp 


T(tj) 


Saoitj), 


S(tj)<Sao(tj) 
S(t  j)>Sa0(t  j) 


(19) 

To  predict  the  amount  of  mine  burial,  BD,  it  is  assumed 
that  the  mine  immediately  settles  to  the  bottom  of  the  scour 
pit  as  the  pit  grows.  Thus,  Bp)(t j)  is  the  maximum  value  of 

S(t<tj ),  which  will  be  denoted  max  (.S’ (f  <  tj  ))•  if  no 

infilling  occurs.  The  depth  of  the  scour  hole,  however, 
depends  on  infill  once  the  currents  subside.  If  the  sediment  is 
assumed  to  fill  the  hole  as  soon  as  the  currents  subside,  the 
function  Sit,)  gives  the  instantaneous  scour  depth  at  time  t  = 
tj.  When  0<S(tj)<Bp),  then  the  scour  process  is  re¬ 
exposing  the  mine.  Thus,  the  burial  of  the  mine  under  this 
assumption  is  given  by 


Bp)(t j)  =  max(S(t j  <f)]-«S(f j)  (20) 


A.  Experimental  Set-Up 

The  mine  was  designed  by  NRL  and  OMNI  technologies 
of  New  Orleans  based  on  an  instrumented  mine  developed  by 
Ingo  Stender  of  Forschungsanstalt  der  Bundeswehr  fur 
Wasserschall-  und  Geophysik  (FWG)  in  Kiel,  Germany  [12], 
Heading  (+/-  1°)  is  measured  with  three  solid-state  compasses 
and  roll  and  pitch  (+/-  1°)  of  the  mine  are  measured  with  a 
three-axis  accelerometer.  Burial  is  measured  by  three  rings  of 
paired  optical  sensors  externally  mounted  at  15°  intervals 
around  the  mine.  Transmitting  optical  sensors  are  LED’s  and 
receiving  sensors  are  phototransistors.  Burial  is  measured  by 
blockage  between  these  72  pairs  of  sensors.  The  mine  is 
made  of  aluminum  and  is  1.5  m  long  and  0.47  m  in  diameter. 
The  weight  in  air  is  619  kg  and  in  water  is  357  kg,  but  these 
weights  are  adjustable.  For  the  Scripps  Pier  experiment,  data 
was  stored  on  a  hard  drive  inside  the  mine.  During  the 
MVCO  experiments,  measurements  were  made  every  2-5 
minutes  and  transmitted  to  a  shore-side  computer. 


Fig.  1.  The  NRL  instrumented  mine  deployed  at  Martha’s 
Vineyard  Coastal  Observatory 
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B.  Numerical  Techniques 

The  scour  model  was  coded  to  run  on  a  PC  in  MATLAB. 
The  code  used  was  adapted  from  the  original  version  written 
by  Friedrichs  [3,  11].  Time  series  inputs  required  for  the 
models  are  orbital  velocity  and  wave  period;  currents  are 
assumed  to  cause  an  insignificant  effect.  For  the  Scripps  Pier 
deployment,  orbital  velocity  was  calculated  from  wave  height 
and  period  data  using  linear  wave  theory  [8].  The  MVCO 
model  runs  used  velocity  data  and  wave  period  measured  by 
ADCP  and  ADV  sensors.  In  addition,  a  sector-scan  sonar 
(SSS)  was  deployed  at  MVCO  for  continuous  recording  of 
acoustic  imagery  of  the  mine  and  surrounding  seabed  [13], 

The  time  step  was  chosen  to  be  20  minutes.  If  data  was 
missing  at  appropriate  intervals,  spline  interpolation  was  used 
to  fill  in  missing  values.  For  the  MVCO  runs,  two  sets  of 
calculations  were  made.  One  used  wave  period  data  to 
calculate  the  wave  friction  factor  given  by  (13).  The  other  run 
used  a  constant  wave  friction  factor  of  0.008.  The  Scripps 
Pier  model  run  used  the  calculated  wave  friction  factor  only. 
In  all  cases,  burial  was  assumed  to  be  no  less  than  10%. 

C.  Scripps  Results 

The  instrumented  mine  was  deployed  off  the  Scripps 
research  pier  in  8  meters  of  water  from  25  July  to  19 
September,  1999  [14],  Wave  energy,  although  low,  was 
sufficient  to  induce  scour.  Figure  2  (top  part)  shows  the  time 
series  data  for  mean  wave  height  and  the  deduced  mean 
orbital  velocity.  The  bottom  shows  the  burial  data  obtained 
by  the  average  blockage  of  light  from  the  3  sets  of  rings. 

Both  prediction  and  data  show  that  burial  occurs  just 
after  a  peak  in  the  orbital  velocity.  The  theory  assumes  that 
the  pit  immediately  refills  after  peaks  in  the  orbital  velocities, 
whereas  data  suggest  that  burial  occurs  at  a  much  slower  rate, 
often  over  a  period  of  3-10  days.  Nevertheless,  the  overall 
level  of  burial  appears  to  roughly  match  the  experimental 
data. 


Scnpps  Pier  NRL  Deployment.  1999 
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Fig.  2.  Mine  burial  data  and  model  output  from  the 
Scripps  deployment,  Summer  1999.  Top:  Oceanographic 
data.  Bottom:  mine  burial  data  and  model  output. 


A  few  issues  became  apparent  during  this  deployment. 
The  model  assumes  that  the  mine  scours  equally  about  the 
surface  and  will  sink  levelly  into  the  sand.  Bed  armoring, 
however,  was  evident  under  the  end  at  ring  3  resulting  from 
winnowing  of  sand  and  scavenging  of  cobble-sized  gravel. 
This  bed  armoring  made  further  scour  improbable  on  that  end 
and  skewed  the  overall  results  toward  less  burial.  In  addition, 
onshore  migration  of  sand  after  the  first  week  increased 
burial.  The  model  does  not  account  for  these  processes. 

I).  Martha 's  Vineyard  Result 

Two  mine  deployments  were  performed  at  MVCO, 
which  is  located  offshore  of  the  southern  coastline  of 
Martha’s  Vineyard.  The  depth  at  the  mine  was  12  meters.  The 
first  deployment  (Fig.  3)  took  place  from  5  December  2001  to 
22  January  2002  [13].  The  second  (Fig.  4)  took  place  from  5 
April  to  28  May  2002. 

Mine  burial  predictions  for  each  deployment  follow 
experimental  results  reasonably;  an  exception  is  notable  at  the 
end  of  the  first  MVCO  deployment  where  more  re-exposure 
is  predicted  than  is  observed.  From  an  examination  of  the 
SSS  imagery,  this  re-exposure  appears  to  have  been  caused 
by  10  cm  changes  in  bed  elevation. 

The  first  deployment  was  the  more  energetic  of  the  two, 
with  complete  burial  achieved  during  a  stormy  event  that 
occurred  between  7  January  and  11  January.  Complete  burial 
was  confirmed  by  the  SSS  and  diver  observations.  The  mine 
only  buried  to  60%  during  the  second  deployment  but 
problems  with  interpreting  the  data  from  ring  3  may  have 
compromised  the  results  and  the  burial  may  be  more  than 
reported  in  Fig.  4.  Unlike  the  Scripps  data,  however,  mine 
burial  occurs  much  more  rapidly,  agreeing  better  in  this 
regard  with  the  model  output. 

IV.  DISCUSSION 

These  experiments  show  that  the  HR  Wallingford 
equations  appear  to  be  promising  for  predicting  the  scour  of  a 
free  body  cylinder  in  sand  with  grain  sizes  on  the  order  of 
200  microns.  They  also  show  that  other  burial  processes  can 
significantly  influence  overall  mine  burial.  These  include 
changes  in  beach  slope  (onshore  sediment  transport),  bed 
armoring  at  the  Scripps  pier  and  migrating  sand  ridges  at 
MVCO. 

Additional  experiments  on  coarse  sand  and  on  sediments 
with  mixtures  of  sand  and  mud  are  needed  to  validate  the 
model.  Additional  observations  by  acoustic  methods  (SSS) 
are  needed  to  gain  greater  insight  into  the  development  and 
size  of  scour  pit  and  rates  of  infilling.  In  addition,  work  is 
required  to  assess  the  applicability  to  a  wider  variety  of 
scenarios.  These  deployments  only  tested  the  applicability  of 
a  cylinder  on  fine  sands.  Experiments  on  coarser  sands  also 
need  to  been  run  and  analyzed.  In  addition,  the  geometrically 
determined  parameters,  p  from  (1)  and  A  and  B  from  (16),  are 
given  only  for  5:1  cylinders,  pipelines  and  pilings.  Work  is 
required  to  determine  what  these  parameters  may  be  for  other 
mine  geometries.  The  usefulness  of  including  the  tidal 
currents  also  needs  to  be  addressed. 
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Fig.  3.  Mine  burial  data  and  model  output  for  the  first 
deployment  at  MVCO,  Winter  2001-2002.  Top:  Oceano¬ 
graphic  data.  Bottom:  mine  burial  data  and  model  output. 


Fig.  4.  Mine  burial  data  and  model  output  for  the  second 
deployment  at  MVCO,  Spring  2002.  Top:  Oceanographic 
data.  Bottom:  mine  burial  data  and  model  output. 


V.  SUMMARY 

The  HR  Wallingford  equations  have  been  applied  to 
predicting  cylindrical  mine  burial  in  fine  sands  and  tested 
against  data  from  three  experiments.  In  these  experiments,  the 
predictions  and  data  generally  match.  Problems  have  been 
identified  for  future  work  on  the  model’s  development.  In  its 
present  state,  it  appears  mature  enough  for  use  as  a  prediction 
tool  for  cylinders  in  fine  sands. 

ACKNOWLEDGEMENTS 

We  are  grateful  for  funding  provided  by  the  Office  of 
Naval  Research  and  Naval  Research  Laboratory  base  funding 
(Program  #  62435N).  We  thank  Carl  Friedrichs  of  the 
Virginia  Institute  of  Marine  Science  (College  of  William  and 
Mary)  for  his  creation  and  supply  of  the  computer  code  for 
the  scour  model  and  helpful  discussions,  Kevin  Briggs  of 
NRL-Stennis  for  plotting  and  providing  the  experimental  data 
in  Fig.  4,  and  Scott  Jenkins  of  the  Scripps  Institute  of 
Oceanography  for  assistance  with  the  Scripps  deployment. 

REFERENCES 

[1]  M.D.  Richardson  and  K.B.  Briggs.  Seabed-Structure 
Interactions  in  Coastal  Sediments.  Proceedings  of  the  4th 
International  Symposium  on  Technology  and  the  Mine 
Problem.  Naval  Postgraduate  School,  Monterey, 
California,  13-16  March  2000. 

[2]  M.R.  Carstens  and  C.S.  Martin  “Settlement  of  Cylindri¬ 
cal  Mines  into  the  Seabed  under  Gravity  Waves”  Geor¬ 
gia  Institute  of  Technology,  Final  Report  A-628,  1963. 

[3]  C.T.  Friedrichs,  “Existing  simple  models  for  mine  burial 
in  sand  by  wave  scour”  Oral  Presentation,  ONR  Mine 
Burial  Program  Modeling  Workshop,  St.  Petersburg,  FL, 
February  2002. 


[4]  R.  Soulsby,  Dynamics  of  Marine  Sands:  A  Manual  for 
Practical  Applications,  Thomas  Telford:  London,  1997. 

[5]  R.J.S.  Whitehouse,  Scour  at  Marine  Structures:  A 
Manual  for  Practical  Applications,  Thomas  Telford: 
London,  1998. 

[6]  B.M.  Sumer,  J.  Fredspe  and  N.  Christiansen,  Scour 
around  vertical  pile  in  waves,  Proc.  Am,  Soc.  Cev. 

Engrs.,  J.  Waterway,  Port,  Coastal,  and  Ocean  Engng., 
Vol.  188,  1992,  pp.  15-31. 

[7]  Reference  [4],  p.  160. 

[8]  R.L.  Soulsby  and  R.J.S.  Whitehouse,  Threshold  of 
sediment  motion  in  coastal  environments.  Proc.  Pacific 
Coasts  and  Ports  ’97  Conf,  Christchurch,  U.  of 
Canterbury,  New  Zeland,  1997,  pp.  144-154. 

[9]  Reference  [3],  Ch.  4-6. 

[10]  J.  Fredspe,  B.M.  Sumer  and  M.M.  Arnskov,  Time  scale 
for  wave  current  scour  below  pipelines,  Int.  J.  Offshore 
Polar  Engng.,  Vol.  2,  1992,  pp.  13-17. 

[1 1]  C.T.  Friedrichs,  Virginia  Institute  of  Marine  Science, 
private  correspondence. 

[12]  S.  Griffin,  J.  Bradley,  M.D.  Richardson,  K.B.  Briggs  and 
P.J.  Valent.  Instrumented  mines  for  mine  burial  studies. 
Sea  Technology,  Vol.  42  No.  11,  2001,  pp.  21-27. 

[13]  M.D  Richardson  and  P  Traykovski.  Real-time 
observations  of  mine  burial  at  the  Martha’s  Vineyard 
Coastal  Observatory.  11  pps..  Proceedings  of  the  5th 
International  Symposium  on  Technology  and  the  Mine 
Problem.  Naval  Postgraduate  School,  Monterey, 
California,  22-25  May  2002. 

[14]  M.D.  Richardson,  P.J.  Valent,  K.B.  Briggs,  J.  Bradley, 
and  S.  Griffin.  NRL  Mine  Burial  Experiments. 
Proceedings  of  the  Second  Australian-American  Joint 
Conference  on  Technologies  of  Mine  Countermeasures, 
Sydney,  Australia,  21-29  March  2001.  Defence  Science 
and  Technology  Organization. 


2354 


